Rabi type oscillations in damped single 2D-quantum dot 
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Abstract 

We present a quantized model of harmonically confined dot atom with inherent damping in the 
presence of a transverse magnetic field. The model leads to a non hermitian Hamiltonian in real 
coordinate. We have analytically studied the effects that damping has on the Rabi type oscillations 
of the system. The model explains the decoherence of Rabi oscillation in a Josephson Junction. 
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I. INTRODUCTION 



Rabi oscillation [1] is one of the fundamental observations in light matter interaction that 
occurs coherently and nonlinearly Q and which has no classical analogue. The generation of 
coherent superposition of quantum states using ultra short laser pulses and the subsequent 
decoherence due to some inherent damping or interaction with the environment is of great 
interest especially in semiconductor quantum dots due to the prospect of future applications 
3], 0, 0] in quantum information processing and making novel laser devices jfj] . Rabi oscil- 
lations using excitons in single quantum dots [7 
different groups in the past few years 
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10| 1 have been studied successfully by 



14j . Control of the decoherence of Rabi 



oscillation in quantum dot, the mechanism of which is still a matter of investigation, has 
attracted wide attention Q]- There has been rapid progress in experimental control 
of dephasing of coherent states in quantum dots [17|. On the contrary, theoretical studies 
on quantum dots leading to the dynamics in the presence of damping are scarce. Till date 
theoretical studies in dots have been done mainly on the basis of damping that has been in- 
troduced phenomenologically. To the best of our knowledge, no quantum theoretical model 
of the dot has been developed with the inherent damping incorporated in the model. 

We develop a model of a damped quantum dot going beyond the phenomenological de- 
scription used so far. The model, we believe provides some insight into the Rabi dynamics. 
Our analytical results lead to an understanding of the experimental observation of decoher- 
ence in Josepshon Junction [Lsl ]. 



II. MODEL 



We show, in what follows, that the damped one electron dot can be described by the 
eigenstates of a quantum Hamiltonian H that is non hermitian. The artificial atom that we 
have modeled is composed of a single electron confined in 2-D by harmonic potential with 
some inherent damping and a homogeneous magnetic field applied normal to the confinement 
plane. Let us start with the classical equation of motion of the damped harmonic oscillator 
which reads 

m e ~r* + 7~r* + hPf = (1) 



2 



where k is the harmonic force constant and 7 is the damping constant and m e is oscillator 
mass. The system described by equation (1) is known to have a time dependent Lagrangian 



and Hamiltonian 
oscillator 
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2l| . There have been many attempts to quantize the damped linear 



24j but a completely satisfactory solution been elusive. The stumbling 



block has been the lack of a time independent Hamiltonian formalism. Recently, however 
such a formalism has been proposed making a definite progress [25|, [26(. We proposed a 



different strategy that brings a non hermitian Hamiltonian formalism. From equation (1) 
we start by noting that it is immediately possible to write down the Euler-Lagrange equation 
for the dissipative system by defining a velocity dependent force F gen and setting 



d (dL\ _dL 
dt \ dr I dr 



gen 



(2) 



where F gen is defined as the negative derivative of Rayleigh dissipative function with 
respect to r [l9]. 



Fgen — ^(/) 

/ is determined by the damping constant 7 and the velocity (r) as follows: 

f = ~ 2 ir 



(3) 



(4) 



Equations (3) and (4) suggest that the time dependent damping force (F d ) is linearly related 
to the velocity: 

F d = -7 r> (5) 
With equation (3) the Euler Lagrange equation (2) now reads 



d_fdL\ _^ + ^/ _ 
dt\dr) dr dr 

Equation (6) requires that the Lagrangian L is chosen as 

L = -m e r 2 + jrr — -kr 2 



(6) 



(7) 



Clearly the Lagrangian of equation (7) is consistent with the equation of motion of the 
damped harmonic oscillator equation (1). Since the momentum p = %- the modified mo- 
mentum for the damped harmonic oscillator becomes 



p = m e r + 7 r 



(8) 
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Let the damped oscillator have a charge 'q' and let it experience an electric field (E) and a 
transverse magnetic field (B). The Lorentz force acting on it is 



F = q 



E+-(vx B) 



= Q 



m 1 f dA " 



-(v x B) 

c 



(9) 



where E = — V0 ((f) =scalar potential) and B = Vxi (A= Vector potential). 

The electric and magnetic fields bring in additional terms in the Lagrangian (L = L, say) 

where 



1-2 / 

-m e r + "yrr — qcp H — A r 



(10) 



, where q<p = \kr 2 , scalar potential. The modified momentum (p) for the system (described 
byL) 



p = m e 'r f + 7"r* H — A. 



(11) 



The modified momentum p leads to the Hamiltonian (H) of the system represented by a 
single carrier electron in a damped quantum dot as follows: 

1 



H 



2m P 



(m e r +7r H — A J • (m e r + 7r H — A.) 

c c 



(12) 



Taking the cyclotron frequency u c = ^^,the confinement potential qcj) = \m e uj 2 {x 2 + y 2 ) 
and replacing the classical operators by their respective quantum analogues, the quantum 
mechanical Hamiltonian of the system in Cartesian coordinates becomes 



H = 



ti 2 



d 2 



d 2 



2m e \dx 2 dy 2 

„2 



ih'-y ( d d " 
^[l + x— + y— 
m e V ox oy 



d 



d 



^ dx X 'dy 



+ 



^-{x 2 + y 2 ) + ^u; 2 (x 2 + y 2 ) + ^m e u 2 (x 2 + y 2 ). 



(13) 



Transforming from Cartesian to polar coordinates the Hamiltonian changes to 

h 2 



H 



where Q 2 , 



2m e \ dr 2 r dr r 2 d<p 2 



dr 



d(j) 



(14) 



+ 



UJr, 



H is manifestly non-hermitian. H may be thought of as defining a set of eigenstates ip n ,i( r , 0) 
with complex energy E n j if we assume that H obeys the energy eigenvalue equation 



(15) 
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A straight forward series solution of equation (15) (Appendix-A) leads to the quantized 
energy eigenvalues of the damped dot: 

<jJ r l 



E, 



n.l 



+ (2n + l + l)tt - «7(2n + I + 1)Q 



(16) 



where 'n' and '/' are principal, and angular momentum quantum numbers, respectively and 



The energy is clearly complex and the imaginary part of it is related to the dissipating 
energy which is given by 



nl 



-j(2n + l + l)Q 



(17) 



Thus, starting from the classical equation of motion of the damped harmonic oscillator 
quantization has been carried out through a Lagrange-Hamiltonian formalism, where the 
Hamiltonian is non-Hermitian 



27 



28] as expected for a non conservative system 



29 



30]. We 



have described the system in terms of real positional coordinates in contrast with attempts 
to handle the problem in terms of complex coordinate 271 ] . 

Thus, proceeding with the assumption that the system described by the non-Hermitian 
iamiltonian of equation (13) satisfies time-independent Schrodinger equation Hip = Eijj 



27] . we have obtained all the quasi energy eigenstates. 



C n 2 r 2 , |,| 

e 2 r l W 



20F 



(18) 



where L]^ is the Laguerre series and C is the normalization constant. In the absence of 



damping these states merge into Fock-Darwin energy spectrum [311. 132|. while the presence 
of damping makes the energy levels quasi stationary. The important outcome is that for a 
known u and u c comparison of the energy separation between two states as observed from 
experiment and obtained from the expressions with and without damping can lead to the 
realization of the intrinsic damping coefficient of a dot system. For damped dot system the 
energy states are shifted from the energy levels without damping and the shifts are more 
pronounced for stronger damping whereas for greater effective mass of the carrier electron 
the effect of damping is somewhat quenched. Since the non-hermitian Hamiltonian obtained 
for the damped dot has complex eigenvalues that correlate with the energy eigenvalues of 
the dot in the limit of zero damping, it could be interesting to investigate the dynamics of 
the damped dot in response to perturbation by laser light. 



III. DYNAMICS OF DAMPED QUANTUM DOT: 

Let us consider the time-dependent Schrddinger equation for the complex energy eigen 
states of H; 

ih d*nAr,t) = {Er _ lV)n ^ n i ^ t) j (1Q) 

The corresponding wave function is decaying and the probability P(r,t) is proportional 
to \ip n j(r,0)\ e~ ~ . The exponential function accounts for the exponential fall-off of the 
amplitude with time, the first factor being the the amplitude of the initial state which is 
now damped. The intrinsic life time T n j = f n l of these quasi- stationary states are therefore 
determined by the damping coefficient and the quantum number characterizing the states. 

We now consider the two energy levels (g and e) of the damped quantum dot system, the 
two states are designated as ip g and ip e are assumed to be well separated from all other states. 
The system interacts with a laser of frequency lo l and cu a = uj eg is the resonance frequency 
(Fig.l). The effect of perturbation produced by the laser can be treated semiclassically using 
the eigenfunctions of the damped dot ^ Ut i as zeroth order wave function. The perturbed 
Hamiltonian is partitioned into H and V, where the unperturbed dot Hamiltonian H 
is given by the equation (13) and the perturbation in the dipole approximation: V = 
e.rE Q cos(u>Lt)- 

We may now consider the semiclassical perturbation treatment based on the damped wave 
functions of the damped dot already obtained. The time-dependent Schrodinger equation 
for the perturbed system is 

\4f(r,t)) = ~H(r,t)\^f(r,t)) (20) 

while the solution is (k=e,g) 

tt(r,*) = J2C k (t)Mr)e- tuJ * t e~^ t (21) 
k 

Projecting on to the states \e) and \g) and integrating over spatial coordinates in each case 
we arrive at the equations governing the time development of the amplitudes (C g , and C e ); 

iC g = C g {t){tof - iW g ) + Ceit^RV^e-^e-^ (22) 
iC e = C e (t)(u? - iWe) + C fl (t)fi fl ^(t)e^- t e- iw -* (23) 

where the Rabi frequency is defined as = ^-{e\r\g) and the dipole approximation has 
been used. 
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Introducing the transformations C g = [C ff e^+^*] and C e = [C e e^ + >^] 

iC g = C e fi fl l?(t)e<-™» 
iC e = C g n R ^(t)e^ + ^ t 



2(f) 



iujT.ttc—wr.'t 



S detuning frequency 



Invoking the rotating wave approximation we get 



iC„ 



Sl R 



SI, 



C P e-^ lt e iSt 



2 



Equation (26) and (27) can be uncoupled by the standard route, leading to 



C g -( 7 u;i-id)C g + -fC g 








Taking the initial condition that C g (0) = 1 ; C e (0) = we get the solutions 



C n 



S^mt , . A . SlRdt 
cos h i— — sm 



n 



Rd 



e 2 



. SIr Slptdt 
i— — sm 



n 



Rd 



where A = ju^ — i5 and f2 M 
Hence, 



- (7^ - i5f. 



e -7(%i)t e -ift e -*-«t 



. . Sl Rd t 

1— — sm 



n 



Rd 



Hence, the excited state population P e = \C e \ is given by 



Qr 2 . 2 ^Rdt 

sm 



n 2 



(24) 
(25) 



(26) 
(27) 



(28) 
(29) 



(30) 
(31) 



(32) 



(33) 



LRd 2 

The result shows that the contribution of a given state to the evolving wave function (\&) of 
the system at a particular time is given in terms of the damping coefficient and the sum of 
the energies of the two levels coupled by laser light. The coherent temporal oscillations of the 
population in the excited state obtained above matches with the experimental observations 
made by Yu et al 



in Jopsepshon phase qubit. The observed oscillatory behaviour of the 



decaying amplitude reported by them is successfully explained by our model based on the 



description of a damped quantum dot by a non-hermitian Hamiltonian in real Coordinates. 
The probability of being in the state 'k' (e or g) at any given time is therefore given by Pk(t) = 
|Cfc(i)| 2 . For the excited state 'e' Figure 2 shows the nature of the time dependence of P e . As 
expected it is coherently oscillatory and exponentially damped. We note that the equation 
33 was earlier developed by Yu et al 18| as the asymtotic limit of solution of the appropriate 
Lioville equation for the density operator under the rotating wave approximation, and used 
to interpret their experimental observation. We have arrived at the same results based on 
the non-hermitian Hamiltonian. 



IV. CONCLUSION 



In summary the proposed model describes correctly the effects of inherent damping in 
a quantum dot. The amount of dissipating energy in a particular state in a quantum dot 
is naturally related to the damping coefficient. The decoherence of Rabi oscillations shows 
that the rate of decoherence is exponentially related not only to the damping coefficient but 
also to the energy separation between the two levels. Again one interesting point is that the 
inherent life times of all the different states is predictable assuming that the life time of any 
one particular state are known from experiment. We also note that the temporal coherent 
oscillation of population in Josephson junction is correctly explained by the present model. 



V. APPENDIX 



In atomic units the Hamiltonian of equation (14) reads 
H -- 



1 ( d 2 Id Id 2 
+ + 



2m P \ dr 2 r dr r 2 d(b 2 



Z7 / ^ ^ d 
m e \ dr 



iu c ( d " 

~2~ \M, 



2„2 



(34) 



where Vi\ = \m e 



Substituting, /(r, <p) 
Schrodinger equation 



, i 2 i , ,2 



2- 



and multiplying both sides by y/2ne '^rs leads to radial 



1 



2m e 



d 2 1 ! 2 n22 



f(r) + ^ - ^ + rf'(r) - i/(r)} = Ef{r) (35) 



or 



+ (t - l2 )\ ~ n d r2 ~ m e^d + 2i7(J + r-^-) + 2Em £ 



dr 2 4 



dr' 



f(r) = 



(36) 



S 



Where ^ = m^ + i + Uo 2 ) 

Substituting r = — f= the radial function f(r) changes to g(x) 



dx 2 



+ T" 



'I 

2 +X ^ 



+ 2£m e 



g(x) = (37) 



'd "d 

Asymptotic analysis leads to 

= ^o(^)^(x)5( 00 (a;) ; # (z) = e - * 2 / 2 , ^(x) = x^+l'l 
g(x) = e~ x2 / 2 V(x)x^ + ^. 

Where V(x) = Y,j fyx^ satisfies laguure series. Again taking z = x 2 the function V(x) 
changes to function q(z), such that q(z) = Y^k a k zk satisfy the Laguure series. 



r d 2 n 1 d J + 1 (u c l-2E), 

\ Z J-I + ( l + l - Z )--{— + m e ,n, > 



4^ 



+2^(1 + +2^— + l- z )]q(z) 
4 



(38) 



where E is complex (= Er — iT). The Lageurre series satisfies equation (36). Accordingly 
the total wave function reads 

^,(r, 0) = — = e^e-^-r' £ 6^** (39) 

2V7T o 
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Figure 1 




r =ya> 

FIG. 1: Energy level diagram for two-level system showing decay rates for ground and excited 
states T g , T e respectively. The excited and ground state frequency difference is denoted by uj eg . 
The Rabi oscillation Qr is introduced by an external laser frequency lol- 



Figure 2 
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FIG. 2: A plot of total population of the excited states versus real time (in a.u) with m e =1 a.u., 
uj c = 10~ 3 a.u., ujq =0.0141 a.u., Qr = 0.01 a.u., (u g + 0J e ) = 0.04 a.u. and 7 = 0.0001 a.u. 
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